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Abstract. A Monte Carlo model has been developed to study the degra- 
dation of <1000 eV electrons in an atmosphere of C0 2 , which is one of the 
most abundant species in Mars' and Venus' atmospheres. The e-C02 cross 
sections are presented in an assembled set along with their analytical rep- 
resentations. Monte Carlo simulations are carried out at several energies to 
calculate the "yield spectra" , which embodied all the information related to 
electron degradation process and can be used to calculate "yield" (or pop- 
ulation) for any inelastic process. The numerical yield spectra have been fit- 
ted analytically resulting in an analytical yield spectra (AYS). We have cal- 
culated the mean energy per ion pair and efficiencies for various inelastic pro- 
cesses, including the double and dissociative double ionization of C0 2 and 
negative ion formation. The energy distribution of the secondary electrons 
produced per incident electron is also presented at few incident energies. The 
mean energy per ion pair for C0 2 is 37.5 (35.8) eV at 200 (1000) eV, com- 
pared to experimental value 32.7 eV at high energies. Ionization is the dom- 
inant loss process at energies above 50 eV with contribution of ~50%. Among 
the excitation processes, 13.6 eV and 12.4 eV states are the dominant loss 
processes consuming ~28% energy above 200 eV. Around and below ioniza- 
tion threshold, 13.6 eV, 12.4 eV, and 11.1 eV, followed by 8.6 eV and 9.3 eV 
excitation states are important loss processes, while below 10 eV vibrational 
excitation dominates. 



1. Introduction 

Photoionization is the main source of electrons and ions in the dayside upper atmo- 
sphere of planets. Photoelectrons, generated due to photoionization process, can have 
enough kinetic energy to ionize the atmospheric constituents and produce secondary elec- 
trons. Similarly, energetic electrons precipitating along the magnetic field lines into the 
auroral atmosphere of planets can ionize the medium producing secondary electrons. Be- 
sides ionization, the electron energy is lost in excitation, attachment, and dissociation. 
Hence, the study of electron energy deposition in atmosphere is an important aspect in 
understanding processes like aurora, dayglow, nightglow [e.g., Bhardwaj and Gladstone, 
2000; Fox et al, 2008]. To model the electron energy degradation in an atmosphere one 
has to first compile cross sections for various loss processes, and then develop an electron 
energy apportionment method, which will distribute the electron energy among different 
loss channels. 

The study of the electron energy degradation in C0 2 is of fundamental interest in 
various fields of science. C0 2 is one of the most important molecules in our solar system. 
It comprise more than 90% of the atmospheres of Venus and Mars. It is also used in 
lasers, gaseous discharge or low power plasma device. Electron energy degradation in 
CO2 gas has important applications to Mars and Venus. Earlier results from Mariner 
satellites and Mars 3 and Mars 4 spacecrafts have confirmed the presence of an ionosphere 
on Mars, and also detected various emission features on Mars [e.g., Barth et al., 1971; 
Dementyeva et al, 1972], which have been studied in detail by recent SPICAM ultraviolet 
spectrometer observations aboard Mars Express [e.g., Bertaux et al, 2006; Leblanc et al, 
2006]. Emissions from Venus have been studied quite extensively by Pioneer Venus [e.g., 
Fox and Bougher, 1991] and by the ongoing Venus Express [e.g., Bertaux et al, 2007]. 
Electron impact excitation and dissociative excitation of CO2 are the key processes in the 
production of several emissions on Mars and Venus. 

In this paper we present a Monte Carlo model which describes the energy degradation 
of <1000 eV electrons in an atmosphere of C0 2 . Earlier studies of electron degradation 
in C0 2 have been carried out by Sawada et al. [1972], Green et al. [1977], and Fox 
and Dalgarno [1979]. Monte Carlo methods are class of numerical methods based on 
stochastic technique. Though it is time consuming, but due to its probabilistic nature, 
it is an excellent technique for studying the energy degradation of particles, provided 
sufficient sample size is taken. Hence, Monte Carlo methods have been widely used in 
problems dealing with energetic particle degradation in gases and in applications to the 
planetary atmospheres [e.g., Cicerone and Bowhill, 1971; Ashihara, 1978; Green et al, 
1977, 1985; Singhal et al, 1980; Singhal and Green, 1981; Singhal and Bhardwaj, 1991; 
Bhardwaj and Singhal, 1993; Michael and Bhardwaj, 2000; Bhardwaj and Michael, 1999a, 
b; Shematovich et al, 2008]. 

In section 2, we present a compilation of all the e-C02 loss processes cross sections 
available up to the present date and fitted them with a simple analytical form. These 
analytically fitted cross sections can be easily used in the Monte Carlo model, which is 
presented in section 3. The output of the Monte Carlo simulation is employed to generate 
a "yield spectrum," which is presented in section 4. The concept of the yield spectrum was 
first introduced by Green et al. [1977] and further developed by many workers [e.g., Green 
and Singhal, 1979; Singhal and Green, 1981; Singhal and Haider, 1984; Green et al, 1985; 



Singhal and Bhardwaj, 1991; Bhardwaj and Singhal, 1993; Bhardwaj and Michael, 1999a]. 
The yield spectra embodied the information about the electron degradation processes and 
can be used to calculate "yield" for any inelastic event. The numerical yield spectrum is 
represented in an analytical form resulting in an analytical yield spectrum (AYS). The 
AYS and its comparison with the numerical yield spectrum is also presented in the section 
4. In sections 5 and 7, we present the calculated mean energy per ion pair and efficiencies 
for inelastic processes, respectively, using AYS and compare them with that obtained by 
using numerical yield spectra. The energy distribution of secondary and tertiary electrons 
produced during ionization events is presented in section 6. Summary of the paper is 
presented in section 8. 

2. Cross sections 

2.1. Total 

The laboratory measured total scattering cross section (TCS) is available between 0.1 
eV and 5000 eV. The TCS for e-CC>2 collision has been measured by several authors in 
different energy ranges - Ferch et al. [1981] in the energy range 0.007-4.5 eV, Buckman et 
al. [1987] 0.1-5 eV, Szmytkowski et al. [1987] 0.5-3000 eV , Kimura et al. [1997] 0.8-500 
eV, Kwan et al. [1983] 1-500 eV, and Garcia and Manero [1996] 400-5000 eV. At low 
energies, the TCS of Szmytkowski et al. [1987], Buckman et al. [1987], and Ferch et al. 
[1981] are in agreement to within 10%. Recently, Zecca et al. [2002] have determined the 
best value of TCS. In the lowest energy range (<1 eV) Zecca et al. [2002] adopted the 
experimental data of Ferch et al. [1981] and Buckman et al. [1987], which are in good 
agreement with each other. In the 1-1000 eV energy range, Zecca et al. [2002] averaged 
the cross sections obtained by Szmytkowski et al. [1987], Kimura et al., [1997] and Kwan 
et al. [1983], with equal weight, to obtain the recommended values, which are in good 
agreement with Garcia and Manero [1996] at higher (>400 eV) energies. In his review, 
Itikawa [2002] has recommended the TCS of Zecca et al. [2002]. The TCS reaches a 
maximum value of 60 x 10~ 16 cm 2 at 0.1 eV [Ferch et al., 1981; Buckman et al., 1987], it 
then goes through a minimum of 5.5 x 10~ 16 cm 2 at 1.9 eV [Szmytkowski et al. 1987]. At 
lower energies a resonant structure is present ~3.8 eV. 

2.2. Elastic 

2.2.1. Differential elastic 

The differential elastic scattering cross section (DCS) for e-C0 2 collision has been mea- 
sured by many authors [cf. review by Itikawa, 2002; Karwasz et al, 2001]. 

In the 1-4 eV energy, the DCS values of Gibson et al. [1999] and Tanaka et al. [1998] 
are in good agreement at forward angles (<50°), however at larger angles they differ by 
20-30%. Overall, at most of the energies there are good agreement in shape between these 
two DCS. At 30, 40, and 50 eV, the DCS measurements of Gibson et al. [1999], Kanik 
et al. [1989], and Tanaka et al. [1998] are in reasonable accord, within the uncertainties 
of each measurement, and at 50 eV the DCS of Gibson et al. [1999] and Register et al. 
[1980] are consistent. At 100 eV, the measured DCS values of Iga et al. [1999] are in good 
agreement with Kanik et al. [1989] and Tanaka et al. [1998]. 

We have taken the DCS values from Tanaka et al. [1998] in the 1-100 eV range, however 
values at 40, 50, 70, 80, and 90 eV are taken from Kanik et al. [1989] which agree well 
with the cross section of Tanaka et al. [1998] in the entire energy range. The DCS values 



in 200-400 eV range are taken from Iga et al. [1999], and those in 500-1000 eV are taken 
from Iga et al. [1984]. In Table 1, we present the DCS values used in this work. 
2.2.2. Total elastic 

Based on the DCS measured by Register et al. [1980], Tanaka et al. [1998] and Gibson 
et al. [1999], Buckman et al. [2002] have determined the total elastic cross section in 
1-100 eV range with an estimated uncertainty of ±30%. Shirai et al. [2001] have reported 
the recommended elastic cross section up to 1000 eV by considering the beam data of Iga 
et al. [1999]. Itikawa [2002] has recommended the elastic cross section of Buckman et al. 
[2002] in the energy range 1-60 eV, and Shirai et al. [2001] in the energy range 100-1000 
eV. The two data sets merge smoothly . 

We have taken total elastic cross section as recommended by Itikawa [2002] . The total 
elastic cross section is fitted using the semi-empirical formula [Bhardwaj and Michael, 
1999a]: 
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where A 1 ,Bi,A 2 , and B 2 are the fitting parameters, whose values are 8.090 x 10~ 16 A , 
2.184 x 10~ 2 A -2 keV, 0.92 A -2 and 5.0 x 10~ 4 A -2 keV, respectively, and E is the energy 
of the electron in eV. Lower limit of fit is 30 eV, and fitted cross section is shown in Figure 
1. At energies below 30 eV it is difficult to fit the cross section using above equation due 
to resonance structure present at low energies (~4 eV), and hence these values are fed 
numerically in the Monte Carlo model. 

2.3. Dissociative electron attachment 

The dissociative attachment process in e-C02 collisions, which mainly occurs at energies 
<12 eV, leads to the formation of negative ions O - , O2 , and C~. Rapp and Briglia [1965] 
measured absolute values of the total cross section for the production of negative ions 
from C0 2 . Orient and Srivastava [1983] obtained the cross section for the production of 
ions and showed that it is the dominant anion. Their values are in agreement with 
those of Rapp and Briglia [1965] within the uncertainty of the cross sections (±20%) and 
the energy scale (±0.1 eV). Spence and Schulz [1974] measured the cross sections for the 
production of C~ and O2 ions. The cross section for production has two peaks of the 
order of 10~ 24 cm 2 at 11.3 and 12.9 eV, while cross section for C~ production has three 
peaks with the largest value of ~2 x 10~ 21 cm 2 . The cross sections for and C~ are 
small compared to that of O - , and hence are not considered in our study. 

We have adopted the cross section values of Rapp and Briglia [1965] for the production 
of 0~ ions from C0 2 . The cross section shows a double-peak structure - peaks at 4.1 
and 8.3 eV, with the later peak value (4.28 x 10~ 19 cm 2 ) about 2.5 times the value of the 
former peak. The cross section for each peak has been fitted with the following analytical 
form [Bhardwaj and Michael, 1999a]: 



rr(F\ ^l 11 (9) 



Here t = (E — W p )/U, where W p is the energy at the peak. The values of the overall 
normalization parameter A and the effective width parameter U for each of the peaks 



along with the parameter W p and threshold energy Wth are presented in Table 2. The 
fitted cross sections along with laboratory measurements are given in Figure 2. 

2.4. Ionization 

The ionization and dissociative ionization of CO2 by electron impact produce singly 
and doubly ionized ions (COj, CO+, C + , + , C++, ++ , and CC# + ). The cross sections 
for these processes have been reported by Rapp and Englander- Golden [1965], Shyn and 
Sharp [1979], Orient and Srivastava [1987], Tian and Vidal [1998], and Straub et al. 
[1996]. Recently, McConkey et al. [2008] have reviewed the electron impact dissociation 
cross sections for C0 2 . For the total ionization cross section, measurements of Orient and 
Srivastava [1987], Tian and Vidal [1998], and Straub et al. [1996] are within the error 
limits with values of Rapp and Englander- Golden [1965] upto 1000 eV, and with the data 
of Shyn and Sharp [1979] in the energy range 50-400 eV. Tian and Vidal [1998] have also 
measured the cross sections for double and triple ionization of CO2 due to electron impact. 
After a survey of the available experimental data, Lindsay and Mangan [2002] suggested 
recommended values of ionization cross section. Their partial cross sections are based 
on measurement of Straub et al. [1996]. For total ionization cross section below 25 eV, 
Lindsay and Mangan [2002] adopted the values of Rapp and Englander- Golden [1965]. At 
energies above 25 eV, they reported uncertainties of 5% for the partial cross sections for 
the production of CO^ , CO + , C + , + , and the total ionization cross section. The cross 
sections at energies below 25 eV have uncertainties of 7%. There are also uncertainties 
in appearance energies of fragmented ions CO + , C + , + , C ++ , and ++ . We have taken 
the appearance energies for the fragmented ions from Itikawa [2002]. 

We have used the dissociative and direct ionization cross sections recommended by 
Lindsay and Mangan [2002] [cf. Itikawa, 2002; McConkey et al., 2008]. The CO^ ion can 
be produced in four excited states, viz., X 2 I1 9 , A 2 n u , B 2 X+, and C 2 SjJ". Cross sections for 
A 2 I1 U and B 2 £+ states have been taken from Itikawa [2002], while the cross sections for 
X 2 I1 3 and C 2 S+ states have been taken from Jackman et al. [1977]. For double ionization, 
cross sections of (CO + ,0 + ), (C + ,0 + ), and (0 + ,0 + ) production have been taken from 
Tian and Vidal [1998] up to 600 eV; these cross sections have not been added in the total 
ionization cross section because they are already accounted in the cross sections for the 
formation of CO + , C + , and + ions. All these cross sections have been fitted using the 
analytical expression [Jackman et al, 1977; Bhardwaj and Michael, 1999a]. 
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Here E is the incident energy in eV, I is the fitting ionization potential in the eV, which is 
generally close to the threshold potential (W t h), and a is in units of 10~ 16 cm 2 . This form 



gives the asymptotic behavior <j(E) oc E~ l lnE at high energies, which is expected from 
the Born approximation. The fitting parameters are presented in the Table 3. The fitted 
cross sections for single and double ionization are shown in Figures 3 and 4, respectively. 

2.5. Excitation cross sections 

2.5.1. Vibrational excitation 

CO2 is a linear triatomic molecule, which has three normal modes of vibration, i.e., a 
bending mode (0 n 0), a symmetric stretching mode (n 0), and an asymmetric stretch- 
ing mode (0 n), with excitation energy 83 meV, 172 meV, and 291 meV, respectively 
[Kochem et al, 1985]. Infrared active (010) bending and (001) asymmetric stretching 
modes in the near-to-threshold region follow the Born approximation. Moreover, the 
structure near the threshold of vibration excitation in CO2 has been investigated by 
Kochem et al [1985], vibrationally inelastic DCS above 4 eV impact energies have been 
measured by Register et al. [1980] for scattering angles 10° — 140° and impact energies 
of 4, 10, 20, and 50 eV, and by Johnstone et al. [1995] for only one scattering angle 
(20°) in the energy region 1 to 7.5 eV. Nakamura [1995] determined the vibrational cross 
section using swarm experiment. Kitajima et al. [2001] made measurements of DCS for 
the electron impact excitation of C0 2 for (010), (100), (001), and (020) vibrational modes 
over the scattering angles 20° — 130° and energy range 1.5-30 eV (except at 4 eV where 
the smallest angle was extended upto 10°), and assigned an uncertainty of 30% to their 
measurements. Their DCS is consistent with the results of previous beam-type measure- 
ments. Itikawa [2002] has extrapolated the DCS of Kitajima et al. [2001] to obtain the 
total vibration cross sections for three modes, which are presented in Figure 1. 

In our studies we have taken cross sections for the three fundamental vibrational modes 
(010), (100), and (001) from Itikawa [2002]. There are other modes also but their cross 
sections are small compare to these three fundamental modes. 

2.5.2. Electronic excitation 

There are several features in the optical and electron scattering spectrum of C0 2 in the 
energy loss range between 7 and 11 eV (Herzberg, 1966; Rabalais et al, 1971; Hall et al, 
1973). Except for Rydberg states, there is still no definite consensus about structure and 
assignment of the excited electronic states of CO2. In the energy loss spectra of CO2, 
Green et al. [2002] have found four clearly distinct peaks at 10.98, 11.05, 11.16, and 
11.40 eV, with an uncertainty of 30% in their results. Itikawa [2002] in his review paper 
has recommended the DCS of Green et al. [2002], for the excitation of the 10.8-11.5 eV 
energy loss states. Recently, Kawahara et al. [2008] have given the integral cross section 
for electronic states X E+ and 1 II U of CO2, based on the DCS measurement of Green et al. 
[2002] in the energy range 20-200 eV. 

Theoretical calculations of electronic structure have also been made by several authors 
[Nakatsuji, 1983; Spielfiedel et al, 1992; Buenker et al., 2000; Lee et al, 1999]. Using 
distorted-wave method, Lee and McKoy [1983] calculated the cross section for the excita- 
tion of eight low lying-states. But there is not much agreement among these calculations. 
In summary, there is still a need for a detailed study of excitation of electronic states of 
CO2 by electron impact. 

We have taken the empirical cross sections of Jackman et al. [1977] for the electronic 
states of C0 2 . These cross sections have been obtained using equation: 
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where qo = 47rao-R 2 and has the value 6.512 x 10" 14 eV 2 cm 2 . The fitting parameter are 
given in Table 4. The parameters for the two states 12.4 and 13.6 eV, that corresponds 
to Cameron band of CO [cf. Sawada et al, 1972] have been modified. The peak cross 
section of their sum is 2.40 x KT 16 cm 2 at 80 eV [Erdman and Zipf, 1983]. 

2.6. Emission 

Electron impact dissociation and ionization of CO2 can result in the production of 
excited fragments of CO, O, and CO2 in the neutral and ionized states, resulting in the 
emissions in the ultraviolet region. These emissions are important for understanding 
phenomena like aurora, dayglow that occur in the atmospheres of Mars, Venus, and CO2- 
containing atmospheres. The strong band systems observed on Mars are Fox-Duffendack- 
Barker bands {A 2 U U ->■ X 2 U g ) and ultraviolet doublet (B 2 £+ ->■ X 2 U g ) of CO^, and 
Cameron bands (a 3 n ->■ X X S+) of CO [Ajello, 1971; Barth et al, 1971; Bertaux et al, 
2006; Leblanc et al, 2006]. Ajello [1971] measured the emission cross sections for the 
A 2 U U ->■ X 2 Ii g and 5 2 S+ ->■ X 2 U g bands of CO^ from threshold to 300 eV. He also 
measured cross sections for the excitation of the fourth positive system of CO (A 1 !! — > 
X 1 S + ), the first negative system of CO + (£? 2 E + — > X 2 S + ) and several atomic multipletes 
of carbon and oxygen produced from dissociative excitation of CO2. 
2.6.1. Emission from CO^" 

McConkey et al. [1968], Ajello [1971], and Tsurubuchi and Iwai [1974] have detected 
emissions corresponding to the following transitions: 



The peak value of cross sections measured by the three groups for the above transitions 
are in good agreement with each other. These emissions are well known in the Mars 
upper atmosphere. Both the ground and excited states of CO^ are known to be linear 
[Herzberg, 1966]. The cross section of Ajello [1971] has too steep an energy dependence 
near threshold compared to McConkey et al. [1968] and Tsurubuchi and Iwai [1974]. In 
his review, Itikawa [2002] recommended the cross sections of Tsurubuchi and Iwai [1974], 
for which the peak values are (8.0 ± 2.0) x 10~ 17 cm 2 at 160 eV for the A — X transition, 
and (4.7 ± 1.2) x 10 -17 cm 2 for the B — X transition. We have taken the cross sections 
for A — X and B — X emissions of CO^ from Itikawa [2002] . These cross sections have 
been fitted using equation (3). The fitting parameters are given in Table 3, and fitted 
cross sections in Figure 3. 
2.6.2. Emission from CO+ 

Only Ajello [1971] has measured the cross section for the emission of first negative 
system (£? 2 £ + — > X 2 S + ) of CO + . The cross section exhibits an appearance potential of 
25.11 eV, and the peak value of cross section is 1.9 x 10~ 18 cm 2 around 100 eV. The cross 
section for the excitation of the first negative system of CO + from electron impact on CO2 
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at 218.9 - 226.8 nm 



is about a factor of 25 less than for excitation of the same system from CO [Ajello, 1971]. 
We have adopted the cross section of Ajello [1971], which has been fitted analytically 
using equation (4); the fitting parameters are given in Table 4. Figure 3 shows the fitted 
cross section along with experimental cross section. 

2.6.3. Emission from CO 

Cross sections for the production of Cameron band system (a 3 II — > X 1 S + ) and fourth 
positive system (A 1 !! — > X X S + ) of CO have been measured by Ajello [1971]. The emission 
cross section for the fourth positive system is very weak and Ajello could not measure 
the cross section near threshold (13.48 eV). For the Cameron band system, Ajello [1971] 
reported relative magnitudes of the cross section for the (0, 1) band at 215.8 nm. The 
upper state (a 3 Il) of Cameron emission is metastable and has a long radiative lifetime 
(~3 ms) [Gilijamse et al, 2007], and kinetic energies of the CO(a 3 n) fragments are in the 
range of 0-1.2 eV [Freund, 1971]. Erdman and Zipf [1983] measured the total cross section 
for CO (a 3 n — > X X S + ) electronic transition. They estimated the absolute magnitude of 
total Cameron band emission cross section of 2.4 x 10~ 16 cm 2 at 80 eV. The Cameron 
band is the brightest emission feature in the UV dayglows of both Mars and Venus as well 
as an important emission in C02-containing atmospheres, e.g. comets. 

2.6.4. Emission from O and C 

Both Ajello [1971] and Mumma et al. [1972] have reported cross section for the emission 
of the O 130.4 nm triplet from electron impact on CO2, but the measurements are not 
consistent with each other. There are many other atomic emissions produced in e-C0 2 
collisions, but they have very small cross sections [cf. van der Burgt et al, 1989]. Kanik 
et al. [1993] have reported the emission cross sections for O, + , C, C + , CO, and CO + in 
the wavelength region 40 - 125 nm. All the cross sections of Kanik et al. [1993] are less 
than 10~ 18 cm 2 . We have adopted the O I and C I production cross sections of Jackman 
et al. [1977]. 

3. Monte Carlo Model 

The transport of radiation is a natural stochastic process that is amenable to the Monte 
Carlo method due to its probabilistic nature. In the Monte Carlo simulation, modeling of 
an inherently stochastic system is carried out by artificial random sampling. In the present 
work we have developed a Monte Carlo model to simulate the local degradation of 1-1000 
eV electrons in an atmosphere of CO2 gas. The energy bin size is taken as 1 eV throughout 
the energy range. In the simulation we have considered elastic scattering between electrons 
and neutral CO2 molecules, and various inelastic processes like ionization, excitation, 
attachment, dissociation, etc; the cross sections for these processes are described in section 
2. Figure 5 illustrates how an individual electron is treated in the Monte Carlo simulation. 

The initial energy E of the electron is fixed at the beginning of the simulation and 
the direction of movement of the electron (6, 0) is decided with the help of two random 
numbers R\ and R2 [random numbers are uniformly distributed in the range (0, 1)] as 



= cos" 1 (l-2 J Ri), 



(5) 



= 2txR 2 . 



(6) 



The distance to next collision is calculated from 



S = -log(l-R 3 )/na T , (7) 

where R 3 is a random number, n is the number density of the neutral target species (taken 
as 1 x 10 10 cm -3 ), and a? is the total (elastic + inelastic) electron impact collision cross 
section. After generating a new random number R4, the probability of elastic collision 
Pel = OeijoT is calculated, if P e i > R4, elastic collision takes place, if P e \ < R4, the 
inelastic event takes place, and in this case we further test for the type of inelastic event 
that has taken place with the help of another random number. 
For elastic scattering the energy loss is calculated as 
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Here S is the scattering angle in the laboratory frame, v and m are the velocity and 
mass, respectively, of the electron, and M is the mass of the target particle. Differential 
elastic cross sections (discussed in section 2.2.1) are used to obtain the scattering angle 8. 
Differential cross sections are fed numerically in the Monte Carlo model at 28 unequally 
spaced energy points (1.5, 2, 3, 3.8, 4, 5, 6, 6.5, 7, 8, 9, 10, 15, 20, 30, 40, 50, 60, 70, 80, 
90, 100, 200, 300, 400, 500, 800, and 1000 eV) and at 20 scattering angles (0°, 5°, 10°, 
15°, 20°, 30°, 40°, 50°, 60°, 70°, 80°, 90°, 100°, 110°, 120°, 130°, 135°, 150°, 165°, and 
180°). At intermediate energies and angular points the values are obtained through linear 
interpolation. The energy AE is subtracted from the energy of the test particle. After 
the collision, the deflection angle relative to the direction (6>, 0) is obtained as 

cos 6" = cos 9 cos 9' — sin 9 sin 9' cos 0', 



cos <p = (cos 9 cos sin 9 sin <p — sin <fi sin 9 sin + sin 9 cos <fi cos 9)/ sin 9 , 



(9) 



sin = (cos 9 cos sin 9 cos — cos sin 9 sin + sin 9 sin cos 9 )/ sin 9 . 

Here 9' , 4> are the scattering angles. 

In the case of an inelastic collision, the next step is to find whether the event is ionization 
or any of the other type of inelastic collision. If the collision is an ionization event, a 
secondary electron is produced. The energy of the secondary electron T is calculated with 
the help of a random number R as [Bhardwaj and Michael, 1999a] 
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Here ^ is the energy of the incident primary electron before the ionization event. r$, 
T A , T A) T B , and T s are the fitting parameters, and I is the ionization threshold. The 
values of these parameters are given in Table 3. If the energy of secondary electron, 
produced in the ionization event, is more than the lowest cutoff energy (which is 1 eV 
in our simulation) then it is also tracked in a same manner as the primary electron (cf. 
Figure 5). The secondary electrons can also cause ionization, producing tertiary electrons, 
which are treated in a similar way as secondary electrons. In the Monte Carlo simulation 
we also follow tertiary and subsequent electrons. The number of secondary, tertiary, and 
subsequent electrons produced during the ionization events are stored in the appropriate 
energy bins. After the type of collision event has been decided, the appropriate energy 
is subtracted from the energy of the particle. All the collision events are recorded in 
the appropriate energy bins corresponding to the energy of the electron at the time of 
collision. The history (track view) of a particle with each interaction event is traced until 
the electron energy falls below an assigned cutoff value, which is 1 eV. The sample size in 
the present study is 10 6 particles for each simulation. 

4. Yield Spectra 

When all the sampled electrons have been degraded, we get a two dimensional yield 
spectrum, which is a function of the spectral energy E and incident primary electron 
energy E , defined as [Green et ai, 1977]: 

U(E,E ) = ^, (11) 

where N(E) is the number of inelastic collision events for which the spectral energy of 
the electron is between E and E + AE, where AE is the energy bin width, which is 1 eV 
in our model. This yield spectrum is related to the degradation spectrum or equilibrium 
flux f(E, Eq) of Spencer and Fano [1954] by the equation 



U{E,E )=<j T {E)f(E,E ), (12) 

where ot is the total inelastic collision cross section. 

The analytical yield spectrum U(E,E ) embodies the nonspatial information of the 
degradation process. It represents the equilibrium number of electrons per unit energy at 
an energy E resulting from the local energy degradation of an incident electron of energy 
Eo, and can be used to calculate the yield Jj of any state j at energy E with the help of 
following equation: 

rE 

Jj(E )= U(E,E )P j (E)dE (13) 
Jw th 

where Pj{E) = <jj(E)/<jT{E) is the probability of occurrence of the jth process whose 
threshold potential is W t h- The yield for a particular process obtained by using the above 



equation is used in the following sections to calculate the mean energy per ion pair and 
efficiencies for various loss processes. Except at very low energies, yield spectrum U (E, E ) 
and probability of excitation Pj{E) both vary with E in a much simpler manner than do 
f(E,E ) and a j(E). 

For many application purposes yield spectrum obtained by equation (11) is represented 
in the following form: 

U(E, E ) = U a (E, E Q ) H(E — E — E m ) + 5(E - E). (14) 

Here H is the Heavyside function, with E m being the minimum threshold of the processes 
considered, and S(E — E) is the Dirac delta function which allows for the contribution 
of the source itself. In atmospheric and astrophysical applications it is convenient to 
represent U a (E,Eo) in an analytical form [Green et al, 1977]: 

U a (E, E ) = A 1 Q + A^/^l (15) 

Here £ = -E /1000 and e = E/I (I is equal to lowest ionization threshold). A 1 = 
0.027, A 2 = 1.20, t = 0, r = 0, and s = —0.0536 are the best fit parameters. 

We have also tried two other analytical forms given by Singhal et al. [1980] and Green 
et al. [1985]. The form given by Singhal et al. [1980] is: 

U a (E,E )=C + C lX + C 2X 2 (16) 

Here x = E^/(E + L); where tt = 0.585 and L = 1.0 and E is in keV, C = 0.0185, C x = 
5.98, and C 2 = 210.4 are fitted parameters. The analytical form given by Green et al. 
[1985] is: 

U a (E, E ) = C + C^Ek + K)/[(E - M) 2 + L 2 ]. (17) 

Here E k = i? /1000, and C , Ci, K, M, and L are the fitted parameters which are 
independent of the energy. The values of these constant parameters are C = 0.0299, 
Ci = 430, K = 0.0041 keV, M = 0.31 eV, and L = 1.9 eV. 

In obtaining our analytical fits we did not include values of the yield spectra very 
close to E because in this regime yield spectra contain the rapid oscillations known as 
"Lewis effect" [cf. Douthat, 1975]. These oscillations are channels with a finite number 
of threshold energies, so that there are only certain energies near E which an electron 
can acquire. Obviously, no electron can acquire an energy between E and E — E m , and 
that is why the Heavyside function H is inserted in the first term on the right-hand side 
of equation (14). The numerical yield spectrum represented analytically using equations 
(15), (16), and (17) is the two-dimensional analytical yield spectrum (AYS). In our studies, 
we have used the AYS obtained using equation (15), which is presented in Figure 6 along 
with the numerical yield spectra obtained by using (14). It is clear from Figure 6 that the 
analytical spectra represents quite well the numerical yield spectra above the ionization 
threshold; however, at lower energies (below 15 eV) the AYS departs from the numerical 
yield spectra. Similar behavior is seen in the AYS of Green et al. [1977]. 



To overcome this deficiency we introduce an additional function to modify the lower 
energy part of the AYS: 



U b (E,E ) 



EyA e x /A 1 
(l + e*) 2 ' 



(18) 



Here x — (E — A 2 )/A 1 , and A , Ai, and A 2 are the fitting parameters. The values of 
parameters are A = 10.095, A\ = 5.5, and A 2 = 0.9. Equation (18) only affects the 
lower energy (<15 eV) part of the fit. The final AYS is the sum of equations (15) and 
(18) which is shown in Figure 6 at several incident energies: depicting a better fit at lower 
energies (>5 eV) as well as at higher energies. 

Because of the simplicity of function and cost effective computational advantage, the 
AYS technique has been widely used in different planetary atmospheres for various aero- 
nomical calculations, like steady state electron fluxes and volume production rates for any 
ionization or excitation state; the details of the computational technique are described in 
earlier papers [e.g., Singhal and Haider, 1984; Bhardwaj and Singhal, 1993; Singhal and 
Bhardwaj, 1991; Bhardwaj et al, 1990, 1996; Bhardwaj, 1999, 2003; Bhardwaj and Michael, 
1999a, b; Michael and Bhardwaj, 2000; Haider and Bhardwaj, 2005]. 

5. Mean Energy per Ion Pair 

The mean energy per ion pair, fj,j, is defined as the incident energy Eq divided by the 
number of ion pairs produced. It can be expressed as 



where Jj(E Q ) is the population of the jth ionization process obtained by equation (13). 
The quantity mean energy per ion pair is known to approach a constant value at higher 
energies. 

Figure 7 shows the mean energy per ion pair for the ions C0 2 (including the ground 
and excited states), CO + , + , C + , CO2 + , ++ and C ++ along with the mean energy per 
ion pair for neutral CO2, solid symbol represents the mean energy per ion pair for neutral 
CO2 obtained directly from the Monte Carlo simulation at few energy points. 

Mean energy for all the ions decreases very rapidly above their threshold value, but 
after ~100 eV \x declines slowly and at higher energies it becomes almost constant. The 
values of \i for CO^, CO+, 0+ and C+ at 200 (1000) eV are 53.6 (51.2), 403 (415), 263.1 
(247.8), and 626.7 eV (576.2) eV, respectively. The mean energy per ion pair for neutral 
C0 2 gas is 37.5 (35.8) eV at 200 (1000) eV. Fox and Dalgarno [1979] reported a value 
of 33.1 eV at 200 eV for the /i, while Green et al. [1977] obtained a value of 34.7 eV at 
200 eV from their MDEB method. The measured value of the mean energy per ion pair 
in neutral CO2 is 32.7 at high energies [Klots, 1968]. Mean energy per ion pair for X 2 II 9 , 
A 2 n u , B 2 £+, and C 2 S+ states of CO^ at 200 (1000) eV are 112.3 (118.4), 180.3 (156), 
301.5 (266.4), and 1999 (1222) eV, respectively. 

6. Secondary Electron Distribution 

During the degradation process, every time the electron undergoes an ionization collision 
event, a secondary electron is produced. The energy of the secondary electron produced 



H (E ) = Eo/JjiEo), 



(19) 



is calculated using (10). The maximum energy of the secondary electron produced can be 
(E — I)/2, where E is the energy of the colliding electron and / is the ionization potential. 
As mentioned before, secondary and tertiary electrons are also treated in the same manner 
as the primary electrons in the Monte Carlo model. The energy distribution of secondary 
electrons is presented in Figure 8 at several incident energies showing the number of 
secondary electrons produced per incident primary electron. The energy distributions of 
tertiary and quaternary electrons, which are presented only at E = 1000 eV, are much 
steeper than that of secondary electrons. Each incident electron of E = 1000 eV, at some 
point of its energy degradation process, produces at least one secondary or tertiary or 
quaternary electron, whose energy is <7 eV. 

7. Efficiency 

As the electron collide with the atmospheric particles, they lose their energy and finally 
become thermalized. The energy of the colliding electron is divided among the various 
inelastic loss processes. Efficiency means the fraction of incident energy of the electron 
which is eventually deposited in a particular loss channel after the completion of the entire 
degradation process. The efficiency, r]j(E ), of the jth process at incident energy E can 
be obtained as 



Vj(Eo) = -=P Jj(E ) (20) 

We have calculated the efficiencies for all inelastic collisions using numerical yield spectra 
obtained from equation (14) and the AYS [sum of equations (15) and (18)]. 

Figure 9 presents efficiencies of various single ionization events producing CO+, CO + , 
+ , and C + . The CO+ has the maximum efficiency throughout the energy region due its 
higher ionization cross section. At 1000 eV, ~31% energy of the incident electron goes 
into CO+ formation, while 5.9%, 9.8%, and 5.0% energy goes into the production of CO + , 
+ , and C + , respectively. At higher energies (>100 eV), increase in the efficiencies for all 
ions is small, but near threshold it falls very rapidly. At threshold, efficiencies for CO^, 
CO+, + , and C+ are 5.1%, 1.1%, 0.16% and 0.19%, respectively, while at 200 eV these 
are 29%, 6.0%, 9.2%, and 4.6%, respectively. Efficiencies for CO+(A-X), CO+(B-X), and 
first negative band of CO+(B-X) are also shown in Figure 9. At 200 (1000) eV, 12.2 
(11.6)% of incident electron energy goes in to the emission CO+(A-X), while 9.8 (11.4)% 
and 3.0 (3.3)% goes in to the emissions CO+(B-X) and CO + (B-X), respectively. 

Figure 10 shows the efficiencies for double ionization of C0 2 . At 200 (1000) eV, effi- 
ciencies for CO++, ++ , and C++ are 0.56 (0.67)%, 0.052 (0.12)%, and 0.092 (0.14)%, re- 
spectively. We have also calculated the efficiencies for (CO+,0+), (C+,0+), and (0+,0+), 
based on cross sections of Tian and Vidal [1998], whose values are 2.7 (3.1)%, 1.8 (2.4)%, 
and 0.96 (1.1)% at 200 (1000) eV. It is clear from Figures 9 and 10, that efficiencies 
calculated from the model and those obtained by using AYS are in good agreement. 

Efficiencies for various excitation processes are presented in Figure 11. The 13.6, 12.4, 
and 11.1 eV states dominate the excitation events having efficiencies 16 (15)%, 12 (13)%, 
and 4.7 (4.2)% at 200 (1000) eV, respectively. Efficiencies of various line emissions of 
atomic oxygen and carbon are shown in Figure 12. Efficiencies for O I (1304), O I (1356), 
C I (1279), C I (1329), C I (1561), and C I (1657), are 0.12 (0.13)%, 0.27 (0.28)%, 0.084 



(0.089)%, 0.035 (0.030)%, 0.10 (0.093)%, and 0.19 (0.18)%, respectively, at 200 (1000) eV. 
Overall efficiencies calculated from numerical yield spectra and AYS for various emission 
and excitation events are in good agreement. 

In Figure 13, we present a summary picture of the electron energy distribution in CO2 
for all the loss processes grouped into important loss channels. At higher (>50 eV) 
energies ionization is the dominant loss process with energy consumption of ~50%. At 
lower energies (<15 eV), 11.1, 12.4, 8.6, and 9.3 eV loss channels are more important. At 
energies below 10 eV, vibration becomes the main loss channel. We have also shown the 
efficiency for total attachment process, which produces negative ion O - . The efficiency for 
anion O production peaks around 8 eV with a value of 0.8%, while it is 0.15 (0.13)% at 
200 (1000) eV. The total efficiency for double ionization, which results in the production of 
CO^ + , ++ , and C ++ ions, is also depicted in the figure. The double ionization efficiency 
raises sharply above 40 eV, having value of 0.4 (0.7)% at 100 (200) eV. Around 1000 eV, 
double ionization efficiency is 0.9%, which is higher than that of 8.6 and 9.3 eV excitation 
states. On the other hand, at energies >100 eV efficiency for dissociative ionization is 
higher than that of the 13.6 and 12.4 eV states. 

8. Summary 

In this paper we have presented a Monte Carlo model for <1000 eV electron degradation 
in CO2 gas. All the e-C02 collision cross sections are compiled and fitted analytically. 
The analytical cross sections are presented in figures along with the laboratory measured 
cross sections for direct comparison, and the fitting parameters are provided in tables. 
The output of the Monte Carlo model is used to calculate the numerical "yield spectra" , 
which is represented by an analytical form. This analytical yield spectra (AYS) can be 
used in planetary atmospheres to determine various aeronomical quantities. We have 
modified and improved the AYS presented by Green et al. [1977] and Singhal et al. [1980] 
by adding a term that provides a better analytical representation of yield spectra at lower 
(<15 eV) energies. The yield spectra is employed to compute the mean energy per ion 
pair and efficiency of various inelastic processes. The mean energy per ion pair for CO2 is 
found to be 37.5 (35.8) at 200 (1000) eV. The energy distribution of secondary electrons 
produced per incident electron is presented at few incident energies. 

Efficiency is an effective measure to know what fraction of the incoming particle energy 
goes into a particular loss channel. We have presented efficiencies for various inelastic 
events calculated by using the AYS as well as by using the numerical yield spectra ob- 
tained from the Monte Carlo model. Efficiencies obtained by the two methods are in 
good agreement. In addition to major inelastic processes, efficiencies are presented for 
the formation of negative ions, double and dissociative double ionization of CO2, and total 
vibrational excitation in the (100), (010), and (001) states. Since the AYS do not represent 
well the numerical yield spectra at very low (<5 eV) energies, the yield for vibrational 
excitation and attachment processes calculated by the AYS would be approximate. Ion- 
ization is the dominant loss process at higher energies, above 100 eV ~50% energy goes 
into ionization. At energies around and below ionization threshold excitation processes 
become important, and at energy below 10 eV, vibration is the dominant loss channel 
consuming more than 70% energy. The 13.6 and 12.4 eV loss channels are also important, 
at 1000 eV, around 28% of incident particle energy goes in to these states. A part of these 



states represents the emissions of Cameron band system, which is an important emission 
in atmospheres of Mars and Venus as well as on comets (Bhardwaj and Raghuram, 2009, 
in preparation). 

Efficiencies presented in this paper can be applied to planetary atmospheres by folding 
them with electron production rate and integrating over the energy These results will 
be useful in the modeling of aeronomical processes in atmospheres of Mars, Venus, and 
CCVcontaining atmospheres. 
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Table 1. Elastic 



differential cross section for electron impact on C0 2 (in units of 10 16 cm 2 
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(2.295) 
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(1.589) 
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3536 





8831 





6294 
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5367 


4.0 


(2.007) 
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7.1830 
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4960 
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5880 
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9487 





8375 





6622 


10.0 


(39.19) 


(29.96) 


(20.74) 


11.520 


2 


2977 
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5342 
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9926 





7430 





6260 


15.0 
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8430 
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5777 


20.0 


(13.80) 


(11.77) 


(9.743) 


7.7149 
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2200 
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90.0 
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9400 
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Table 1. Contd. 
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(o 
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Values inside the bracket indicates a linearly extrapolated value. 



Table 2. Parameters for Electron attachment process 





Wp 


A 


U 




o~i 


4.3 


0.0013 x 10- ie 


0.22 


3.4 


0-11 


8.1 


0.0056 x 10~ 16 


0.33 


5.9 



Table 3. Parameters for various ionization processes 







I 


K 




J 








T 5 






COf (Total) 


13.76 


13.76 


9.83 


0.0 


40.59 


1.050 


18.61 


-13.23 


-0.847 


875 


44.52 


co+(x 2 n 9 ) 


13.76 


13.76 


3.480 


0.0 


4.099 


-2.35 


11.11 


-13.26 


-0.847 


1000 


27.52 


co+(A 2 n u ) 


17.8 


17.8 


8.632 


6.0 


86.36 


1.004 


12.00 


-18.80 


-1.996 


550 


10.20 


CO+(B 2 £+) 


18.1 


18.1 


4.632 


5.0 


85.36 


1.004 


12.00 


-18.5 


-0.978 


450 


10.20 


C0 2 +(C 2 £+) 


19.4 


19.4 


0.580 


0.0 


21.19 


1.270 


10.98 


-19.00 


-0.887 


1000 


38.80 


CO+ 


24.76 


24.76 


1.347 


15.00 


6.650 


1.256 


9.556 


-24.0 


1.887 


800 


25.52 


0+ 


24.5 


24.5 


2.399 


40.00 


11.34 


-1.10 


13.42 


-24.0 


-0.587 


0.0 


0.0 


C+ 


29.5 


29.5 


1.659 


0.0 


53.54 


0.625 


10.62 


-29.0 


2.800 


21.8 


44.00 


C++ 


79.94 


79.94 


0.0012 


25.00 


0.100 


0.656 


79.11 


0.0 


0.8473 


0.0 


0.0 


0++ 


94 


94 


0.0055 


10.00 


1.700 


-3.156 


20.11 


0.0 


-0.8473 


0.0 


0.0 


CO++ 


37.6 


44.75 


0.0583 


1.20 


0.250 


0.0 


11.57 


35.26 


1.548 


800 


30.52 


(CO+,0+) 


44.7 


44.7 


0.285 


1.200 


0.550 


0.0 


11.97 


35.26 


1.547 


650 


30.52 


(C+,0+) 


44.7 


44.7 


0.288 


1.200 


0.550 


0.0 


10.40 


45.00 


1.547 


750 


0.0 


(0+,0+) 


44.7 


44.7 


0.158 


1.200 


0.450 


0.0 


5.400 


15.00 


5.547 


750 


0.0 



Table 4. Parameters for various excitation and emission processes 



Excitation states 


vv 




Ct 




P 






n 




r 


i\.r . 


[Vibration, (010)1 1 

L i \ / J 


0.080 


2. 


750 


1. 


000 


0.080 





.750 


0. 


,000060 


0.0 


[Vibration, (100)1 

L ' \ / J 


0.180 


1. 


070 


1 


,000 


0.180 





.750 


0, 


.000031 


0.0 


[Vibration, (001)1 


0.290 


2. 


910 


0. 


500 


0.300 





.810 





.000445 


0.0 


8.6 eV state 


8.600 





,556 


2, 


,000 


8.600 





.936 





,060600 


0.0 


9.3 eV state 


9.300 





,603 


2. 


000 


9.300 





.909 





.064000 


0.0 


11.1 eV state 


7.760 


0. 


,246 


3 


,000 


11.100 


1 


.110 


4 


,420000 


0.0 


[12.4 eV state] 


9.610 


0. 


338 


3. 


000 


12.400 





.830 


6. 


.700000 


0.0 


[13.6 eV state] 

L J 


10.50 


0. 


,625 


3 


,000 


13.600 





.849 


3 


,350000 


0.0 


15.5 eV state 


15.50 


0. 


739 


2. 


,000 


15.500 





.793 


0, 


.139000 


0.750 


16.3 eV state 


12.30 


0. 


.605 


3, 


000 


16.300 





.911 


0, 


.716000 


0.750 


17.0 eV state 


13.00 


0. 


,649 


3 


,000 


17.000 





.878 





,114000 


0.750 


17.8 eV state 


14.00 


0. 


,977 


3. 


000 


17.800 





.725 





.051100 


0.750 


OI (1304) 


20.10 


0. 


,599 


3 


,000 


22.000 


1. 


.000 





,127000 


0.750 


OI (1356) 


16.40 


0. 


.600 


3 


.000 


20.400 





.944 


0, 


.168000 


0.500 


CI (1279) 


15.70 


1 


,000 


3 


.000 


26.200 





.643 





,010400 


0.500 


CI (1329) 


21.80 


1 


.000 


3 


.000 


20.900 


1 


.040 





.020200 


0.500 


CI (1561) 


22.40 


1 


.000 


3 


.000 


24.500 





.982 





.053800 


0.500 


CI (1657) 


21.10 


1 


,000 


3 


.000 


24.100 


0. 


.947 





,872000 


0.500 


[CO + (first negative)] 


18.13 





,656 


2.54 


25.11 





.804 




1.055 


0.0 



Parameters are taken from Jackman et al. [1977], except for the states which are 
inside the square brackets whose parameters have been modified. 
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Figure 1. The e-C0 2 total elastic cross section and vibrational excitation cross sections for 
three modes. For elastic cross section, symbol represents the cross section values of Itikawa 
[2002], and solid curve represents the analytical fit using equation (1). Dashed curve represents 
the vibrational excitation cross sections taken from Itikawa [2002] 
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Figure 2. Dissociative electron attachment cross section for the formation of 0~ ion. Symbol 
represents the values of Itikawa [2002] based on Rapp and Briglia [1965]; solid curve represents 
analytical fit of 0~ cross section using equation(2). I and II denotes the first and second peak. 




Figure 3. Ionization and emission cross sections of CO2. Symbol represents the values 
of Itikawa [2002], and solid curve represents the analytical fits using equation (3) except for 
CO + (B-X) state, which is fitted using equation (4). Note that the cross section for CO+(B 2 E+) 
state has been plotted after multiplying by a factor of 3. 




ENERGY (eV) 

Figure 4. Cross sections for electron impact double ionization of CO2. Symbol represents 
the measured cross section values and solid curve represents the analytical fit using equation 
(3). Cross sections for (CO + ,0 + ), (C + ,0 + ), and (0 + ,0 + ) have been taken from Tian and Vidal 
[1998], and that for CO++, ++ , and C++ from Itikawa [2002]. Cross section for 0++ has been 
plotted after multiplying by a factor of 2. 
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Figure 5. A simplified flow diagram of the Monte Carlo simulation. The diagram shows 
flow upto secondary electron, but tertiary and subsequent electrons are also treated in a similar 
manner in the simulation. 
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Figure 6. The numerical yield spectra from the Monte Carlo model (solid curve) and AYS 
using equation (15) (dotted curve) at incident energies (E ) 30, 100, 200, 500, and 1000 eV. 
Dashed curve represents the improved AYS calculated by summing equations (15) and (18). The 
yield spectra at 100, 200, 500, and 1000 eV are plotted after multiplying by a factor of 10, 80, 
600, and 5000, respectively. 
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Figure 7. The mean energy per ion pair for ions CO^ (COj is the sum of four states X 2 IT S , 
A 2 II U , B 2 £+, and C 2 £+), C0+, 0+ C+ C0£ + , 0++, and C++, and the neutral C0 2 gas (total), 
symbol represents the \i calculated using numerical yield spectra for the CO^ and neutral C0 2 . 
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Figure 8. The energy distribution of secondary electrons at six incident energies (E ): 50, 
100, 200, 500, 700, and 1000 eV. N e represents the number of secondary, tertiary, or quaternary 
electrons produced per incident primary electron. Dotted curve and dashed curve represent 
energy distribution of tertiary and quaternary electrons, repectively at E = 1000 eV. 
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Figure 9. Efficiencies of various ionization and emission processes. Symbol represents the 
efficiency calculated by numerical yield spectra and curves represent the efficiency calculated by 
the AYS. CO^(A 2 n u ) and CO^(B 2 S+) represent FDB (A 2 U U ->■ X 2 U g ) and ultraviolet doblet 
(£? 2 £+ — > X 2 Tlg) emissions, respectively, and CO + (B-X) represents first negative band emission 
of CO+ ion. 
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Figure 10. Efficiencies for double ionization of C0 2 due to electron impact. Symbol represents 
the efficiency calculated by numerical yield spectra and curves represent the efficiency calculated 
by the AYS. 
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Figure 11. Efficiencies of various excited states. Symbol represents the efficiency calculated 
by numerical yield spectra and curves represent efficiency calculated by the AYS. 




Figure 12. Efficiencies of various oxygen and carbon line emissions. Symbol represents the 
efficiency calculated by numerical yield spectra and curves represent the efficiency calculated by 
the AYS. 
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Figure 13. Efficiencies for various important loss channels calculated using numerical yield 
spectra. Dissociative ionization includes the production of + , C + , and CO + ions; double ion- 
ization includes the production of CO^ + , ++ , and C ++ ions; and attachment denotes the 
production of 0~ ion. 



